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Superior thermal conductivity and extremely high mechanical strength in 
polyethylene chains from ab initio calculation 
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The upper limit of the thermal conductivity and the mechanical strength are predicted for the 
polyethylene chain, by performing the ab initio calculation and applying the quantum mechanical 
non-equilibrium Green's function approach. Specially, there are two main findings from our calcula- 
tion: (1). the thermal conductivity can reach a high value of 310 W/K/m in a 100 nm polyethylene 
chain at room temperature; (2). the Young's modulus in the polyethylene chain is as high as 374.5 
GPa, and the polyethylene chain can sustain 32.85% ± 0.05% (ultimate) strain before undergoing 
structural phase transition into gaseous ethylene. 
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Polymer is a useful organic material in many impor- 
tant fields such as solar cell,— where high thermal con- 
ductivity will be helpful to deliver heat promptly and 
prevent the device from being broken down by localized 
heating. However, it is well known that the polymer is 
a poor thermal conductor with thermal conductivity less 
than 1 W/K/m due to a huge number of interfaces inside 
the material. 2 In recent years, there have been increasing 
efforts to enhance the thermal transport probability of 
polymer. A straightforward method is to disperse some 
materials of high thermal conductivity (such as carbon 
nanotube or graphene) into polymer, which forms the 
polymer conrposite. This method can only slightly en- 
hance the thermal conductivity of the polymer, because 
the heat transport is still strongly limited by interfaces in 
the composite. 3 Recently, Liu and Yang found that the 
thermal conductivity of the polymer can be considerably 
increased by mechanical strain^ The enhancement of the 
thermal conductivity is attributed to a better alignment 
of polymer chains under mechanical strain. It implies 
that a more efficient way to enhance the thermal conduc- 
tivity is to produce a crystalline polymer, as confirmed by 
a recent study. 5 The properties of the crystalline polymer 
are closely related to the single polymer chain, and thus 
it is meaningful to investigate a single polymer chain. 
Furthermore, it may be possible to prepare a single poly- 
mer chain in the laboratory, considering the rapid de- 
velopment of the experimental detection technique. 6 As 
polyethylene is the most widely applied polymer, it is 
of practical significance to investigate the properties of 
a single polyethylene chain. Henry and Chen applied 
the molecular dynamics simulation to study the thermal 
transport in the polyethylene chaini^ We are more in- 
terested in an accurate prediction for the upper limit of 
the thermal conductivity in the polyethylene chain, since 
it is almost an impossible task to mimic precisely the 
complex experimental set up in a theoretical theme. 

In this paper, we study the thermal transport and me- 
chanical properties of the polyethylene chain by applying 
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FIG. 1: (Color online) Electronic band structure for the 
polyethylene chain. The energy is with reference to the Fermi 
level. Inset shows the relaxed configuration of a unit cell 
C2H4. 



the ab initio calculation and the non-equilibrium Green's 
function (NEGF) approach. The obtained thermal con- 
ductivity serves as an upper limit of the experimentally 
measured thermal conductivity for this system. For a 
polyethylene chain of 100 nm, the predicted thermal con- 
ductivity is as high as 310 W/K/m at room temperature. 
Our calculation shows that a single polyethylene chain is 
of high strength with Young's modulus Y — 374.5 GPa 
and ultimate strain e c = 32.85% ± 0.05%. 

In the ab initio calculation, we use the SIESTA 
package^ to optimize the structure of the polyethy- 
lene chain. The local density approximation is applied 
to account for the exchange-correlation function with 
Ceperlcy-Aldcr parametrization 10 and the double-^ ba- 
sis set orbital (DZ) is adopted. During the conjugate- 
gradient optimization, the maximum force on each atom 
is smaller than 0.005 eV/A. A mesh cut off of 120 Ry 
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FIG. 2: (Color online) Phonon dispersion and Young's mod- 
ulus of polyethylene, (a). Four low- frequency branches of the 
phonon dispersion, (b). The phonon dispersion in whole fre- 
quency range, (c). The transmission from NEGF approach. 



is used. Periodic boundary condition is applied in the 
growth direction, while free boundary conditions are ap- 
plied in the two transverse directions by introducing suf- 
ficient vacuum space. There are two types of calcula- 
tions in this work. In the calculation of the electron and 
phonon band structures, a unit cell C2H4 is used and an 
8x1x1 Monkhorst-Pack k-grid is chosen for the sampling 
of the one-dimensional Brillouin zone. In the investiga- 
tion of the Young's modulus and the thermal transport, 
a finite polyethylene chain of 20 unit cells is calculated 
and the Gamma point k sampling is adopted. 

The inset of Fig.[T]shows the optimized configuration of 
a C2H4 unit cell in the polyethylene chain. The C-C bond 
length is 1.51 A, while the C-H bond length is 1.12 A. The 
two angles are 6> CCC = 113.4° and 6> HC h = 104.8°. These 
structure parameters are in good agreement with the ex- 
perimental results for the crystalline polyethylene»ii~— 
Fig.[T]is the electronic band structure of the polyethylene 
chain, yielding a large band gap of 9.58 eV, which coin- 
cides with the experimental value of 8.8 eV— and other ab 
initio calculation^. The large band gap in the polyethy- 
lene chain is related to the tetrahedral (sp3) structure 
of the carbon atom in this material, inducing a strong 
binding for all electrons. Fig. [5] (a) displays the phonon 
dispersion in the low-frequency range. The wave vector 
k is along the growth direction and in the unit of 2n/a, 
with a = 2.53 A as the size of the unit cell C2H4. The 
twisting (TW) phonon has nonzero frequency at k = 0. 
This is due to the loss of the rigid rotational invariance 
symmetry^ in the ab initio calculation. For the three 
translational acoustic phonons, the sound velocities are 
9.35 km/s for the longitudinal acoustic (LA) phonon, and 
14.14 km/s and 15.03 km/s for the two transverse acous- 
tic (TA) phonons. It is quite interesting that the sound 
velocities for TA phonons are faster than those for the LA 
phonon. This phenomenon implies that the polyethylene 
chain is more stiff in the transverse directions, owing to 
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FIG. 3: (Color online) Thermal transport properties in the 
polyethylene chain, (a). Thermal conductance versus temper- 
ature, (b). Thermal conductivity (k) calculated from ther- 
mal conductance (<r) by: k, — aL/s, with cross-sectional area 
s = 18.24 A. The length L = 100 nm. 



its zigzag configuration formed by C-C bonds along the 
growth direction. It should be noted that these three 
sound velocities are actually quite large and very close 
to that of the single-walled carbon nanotube, where the 
LA and TA sound velocities are around 20 km/s and 
10 km/s, respectively. 18 The full phonon dispersion is 
demonstrated in Fig. [2] (b) for the polyethylene chain. 
There are four typical regions. Below 400 cm" 1 are the 
acoustic phonons, as zoomed in in panel (a). The fre- 
quency of the C-C vibration is around 1300 cm -1 , which 
is quite close to that in diamond,— due to a similar sp3 
bonding in both systems. The C-H bond is character- 
ized by its vibration frequency around 1500 cm -1 . The 
highest frequency around 3000 cm -1 corresponds to the 
vibration of the H atom, which is also observed in the 
Raman measurement (see eg. Ref. [20). 

Fast sound speeds indicate strong thermal transport 
capability of the polyethylene chain. The thermal trans- 
port can be investigated by different methods. Classi- 
cal results can be obtained from the molecular dynamics 
simulation of the thermal transport, where the phonon- 
phonon scattering dominates the transport property 7,8 
From the theoretical point of view, it is a big challenge 
to provide accurate prediction for the thermal conductiv- 
ity, because the samples in the experiment always possess 
various unpreventable defects. Hence, a more practical 
task is to provide an accurate (quantum) prediction for 
the upper limit of the thermal conductivity. For this 
purpose, we can apply the ballistic NEGF method . 21 i 22 
It is based on quantum mechanics. The phonon-phonon 
scattering is ignored in the ballistic transport region, 
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FIG. 4: (Color online) Mechanical properties of the polyethylene chain, (a). The density of strain energy versuses strain. 
The Young's modulus Y = 374.5 Gpa is obtained by fitting strain energy density to E = 0.5Ye 2 in small strain region 
e G [—0.01, 0.01]. A big jump in the curve reveals the ultimate strain to be e c = 32.85% ± 0.05%. Inset shows two neighboring 
ethylene moleculars after the phase transition of the polyethylene chain at e = e c . (b). Stress-strain relation. 



which is actually quite reasonable for nano-materials. 
The combination of the NEGF and the ab initio calcula- 
tion can provide us an accurate upper limit for the ther- 
mal conductivity^ In the NEGF approach, the thermal 
conductance is calculated by the Landauer formula: 
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du>htL>T[uj] 



dn(uj,T) 



dT 



(1) 



where ft is the Planck's constant. n(u,T) is the Bosc- 
Einstein distribution function. The transmission T[u\ is 
obtained from the Caroli formula: 



7>] = Tr(G r T L G a T R ) 



(2) 



where G r is the retarded Green's function. G a = 
(G r ) is the advanced Green's function and T^/jj is 
the self-energy. These Green's functions can be calcu- 
lated from the force constant matrix from the ab initio 
calculation. 23 Fig. H](c) shows the transmission function 
for the polyethylene chain. The function exhibits some 
regular steps, due to the absence of phonon-phonon scat- 
tering. Fig. [3] (a) shows the temperature dependence for 
the thermal conductance. The thermal conductance (ct) 
does not depend on the length of the system. It can be 
used to get the thermal conductivity (k) of a polyethy- 
lene chain with arbitrary length L: n = aL/s, where 
s is the cross-sectional area of the polyethylene chain. 
Fig. [3] (b) gives the thermal conductivity for a polyethy- 
lene chain of 100 nm. Particularly, the room-temperature 
thermal conductivity is as high as 310 W/K/m. The 
thermal conductivity will increase with increasing length 
of the polyethylene chain, which is quite similar as the 
thermal conductivity of two-dimensional graphene^ - — 
We should stress that this value serves as an upper limit 
for the thermal conductivity of the 100 nm polyethylene 



chain. If the experimental samples are of high quality, 
then the measured thermal conductivity should approach 
310 W/K/m from lower side. 

A fast sound speed in the polyethylene chain also im- 
plies high stiffness of this structure. Indeed, we find large 
value for the axial Young's modulus from Fig. 2] (a). 
In the figure, the density of the strain energy is calcu- 
lated for a finite polyethylene chain of 20 unit cells un- 
der axial uniform strain. In the calculation of the vol- 
ume, the cross-sectional area s = 18.24 A 2 is chosen as 
half of the experimental value for the unit cell of crys- 
talline polyethylene with two polyethylene chains>2£l In 
the small strain region e G [—0.01,0.01], the nonlinear 
effect is neglectable, so the elastic parameter can be ob- 
tained by fitting the strain energy to E(e) — l/2Ye 2 
in this strain region. The obtained Young's modulus 
Y = 374.5 GPa. This value is even higher than the 
Young's modulus of crystalline polyethylene around 300 
GPai^Sr— In larger strain region e G [—0.2,0.2], the en- 
ergy curve is obvious asymmetry under tensile and com- 
pressive strain, which illustrates a strong nonlinear effect 
for the polyethylene chain in the large strain region. This 
nonlinear effect holds such property that the interaction 
is stronger for compressive strain than that of the tensile 
strain, which will result in a positive thermal expansion 
phenomenon during heating. With further increase of the 
tensile strain, the strain energy density increases linearly 
and the polyethylene chain is broken at a critical strain 
e c = 32.85% ± 0.05%, which serves as the ultimate strain 
of the polyethylene chain. The jump of the strain energy 
density at e = e c corresponds to a bond energy of 1.031 
eV for each sp3 C-C bond in the polyethylene chain. The 
discontinuity of the energy jump is a characteristic fea- 
ture for the broken of a single chain. In the structural 
phase transition, the sp3 bonds in the polyethylene are 



broken and new sp2 bonds are formed. According to this 
bond transition, the polyethylene chain changes into 20 
gaseous ethylene moleculars. Two ethylene moleculars 
are shown in the inset of Fig. @] (a) after the phase tran- 
sition of the polyethylene chain. The C-C and C-H bond 
lengths are 1.332 A and 1.103 A, respectively. The angles 
are 6 HC c = 122.2° and HC h = 115.7°. These struc- 
ture parameters agree quite well with the experimental 
datai^i The distance between two neighboring ethylene 
moleculars is about 3.356 A, which falls in a typical inter- 
action range of the van der Waals potential. The phase 
transition phenomenon can also be well captured by the 
stress-strain relation shown in Fig.[4](b), where the stress 
reaches a saturate value before e c and suddenly jumps to 
a very small value at e = e c . 

In conclusion, we have performed ab initio calculation 
to study the phonon, mechanical, and thermal transport 
properties of the polyethylene chain. Our results show 



that the sound speeds are very high in this system: 9.35 
km/s for the LA phonon, and 14.14 km/s and 15.03 km/s 
for the two transverse phonons. Using the NEGF quan- 
tum approach, we provide an upper limit for the ther- 
mal conductivity of the polyethylene chain of any length; 
particularly, the room-temperature thermal conductivity 
can be as high as 310 W/K/m for a 100 nm polyethylene 
chain. In consistent with the high speed, we get a large 
Young's modulus (374.5 GPa), which is even larger than 
that of the crystalline polyethylene. After increasing the 
applied tensile strain, we observe the structural transition 
from a polyethylene chain into gaseous ethylene, yielding 
an ultimate strain as high as e c — 32.85% ±0.05% for the 
polyethylene chain. 
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